Emergent heavy fermion behavior at the Wigner-Mott transition 
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We study charge ordering driven by Coulomb interactions on triangular lattices relevant to the 
Wigner-Mott transition in two dimensions. Dynamical mean-field theory reveals a pinball liquid 
phase whose low energy properties are captured by the periodic Anderson model. Strong quasipar- 
ticle renormalization is predicted at the onset of Kondo antiferromagnetic coupling between localized 
electrons (pins) and itinerant carriers (balls). Non- Fermi liquid behavior can occur in such charge 
ordered systems due to spin-flip scattering of itinerant electrons off the pins in analogy with heavy 
fermion compounds. 
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The interplay between charge ordering (CO) and the 
Mott metal-insulator transition (MIT) is an important 
open problem in correlated materials, whose manifesta- 
tions are common in narrow-band systems such as the 
transition metal oxides [l| (cuprates, manganites, nick- 
elates, cobaltates), low-dimensional organic conductors 
y as well as transition-metal dichalcogenides [3|. On 
one hand, charge ordering phenomena are favored by the 
presence of electronic correlations, as these slow down 
the electronic motion and reduce the tendency to charge 
delocalization. Conversely, the reconstruction of the elec- 
tronic spectrum upon entering a phase with broken sym- 
metry can enhance the effects of correlations, extending 
the realm of Mott physics away from the usual integer fill- 
ings [3, 14411 1 . An illustration of this phenomenon is real- 
ized in the dichalcogenide compound lT-TaS2 |3l.ll2l.ll3|. 
where nominally weak electronic correlations are so effec- 
tively enhanced upon formation of a charge density wave 
that they can trigger a Mott metal-insulator transition. 

A fundamental relation between charge ordering and 
the Mott metal-insulator transition has also been pointed 
out recently in the different context of the two- 
dimensional electron gas (2DEG), leading to the concept 
of Wigner-Mott transition [6|-|9[. Close to the Wigner 
crystal phase induced by the long ranged Coulomb re- 
pulsion at low density, the electrons form an effectively 
narrow-band system at half-filling (one electron per site 
of the Wigner crystal JJij), whose physical properties 
therefore resemble those of a Mott insulator. Corre- 
spondingly, it has been shown that the melting of the 
Wigner crystal shares several common features with the 
Mott transition [i,0,[l3. 

In the present Letter, we explore the interplay of charge 
order and correlations starting from the extended Hub- 
bard model on the triangular lattice, which frustrates 
the insulating CO found in bipartite lattices. Based 
on DMFT calculations at T = we find that the 



nearest-neighbor Coulomb repulsion stabilizes a charge 
ordered metallic phase over a broad range of concentra- 
tions. In this phase, termed pinball liquid and sketched in 
Fig. [U quasi-localized electrons (pins) coexist with itin- 
erant electrons (balls), whose dynamics become heavily 
renormalized for sufficiently large carrier concentration. 
The origin of this renormalization is elucidated via a 
mapping to a periodic Anderson model (PAM): tuning 
the hybridization between the two species reveals an an- 
tiferromagnetic coupling leading to Kondo screening of 
the localized moments. Although our calculations are 
limited to the zero temperature case, the strong analo- 
gies with heavy fermions indicate that non-Fermi liquid 
behavior can arise in correlated CO systems on frustrated 
lattices in a broad temperature range above the Kondo 
breakdown temperature scale. 

Model. The reference model to study the interplay 
between charge ordering and electronic correlations is the 
Extended Hubbard Model (EHM): 



i/ = --i ^ i4^Cj^ + h.c.) + UY^ n.^Uii 



(1) 
which describes electrons moving on a lattice with trans- 
fer integral t and interacting via both on-site ([/) and 
nearest-neighbor Coulomb interactions (F). We set the 
energy scale t — \. Studies of the EHM have mostly 
focused on bipartite lattices [2, Ij, la, LZn9| (the Bethe lat- 
tice, or the square lattice in two-dimensions), as these 
can naturally accommodate for charge ordering at a com- 
mensurate occupation of one charge per two lattice sites, 
corresponding to one quarter band filling. There are im- 
portant cases however where the relevant lattice geom- 
etry is non-bipartite, adding frustration to the already 
rich phenomenology of correlated CO. A notable exam- 
ple is the two-dimensional Wigner crystal itself, where 
charges arrange on a triangular pattern in order to min- 
imize their electrostatic energy. How the analogy be- 



tween the Wigner crystal transition and the Mott MIT 
is modified when one considers the appropriate triangu- 
lar geometry 1J| remains an open question, that was not 
addressed in Refs. [0, 01- A second example is the fam- 
ily of two-dimensional organic conductors 0-ET2 [16J, ll7| , 
where the molecular lattice is triangular and the band fill- 
ing is fixed by charge transfer to one hole per two molec- 
ular sites. In this framework, the EHM has been studied 
intensively to address the phenomenon of charge order- 
ing, and fundamental differences between the square and 
triangular lattice have emerged. Most prominently, for 
isotropic V, on the triangular lattice, geometrical frus- 
tration prevents the two-fold CO characteristic of bipar- 
tite lattices, and a three-fold CO is realized instead (cf. 



Fig. [lb) [SElSfi 



DMFT solution. We solve the model Eq. ([T]) using 
dynamical mean field theory (DMFT). In its simplest 
form, this method maps a single-band Hubbard model 
onto an Anderson quantum impurity model in an ef- 
fective metallic bath that is determined self-consistently. 
Allowing for charge ordering [J] with iV-fold symmetry 
breaking implies solving N interconnected (through hop- 
ping amplitudes) impurity models, one for each sublat- 
tice. Calling G° ^ri^^n) the Green's function of the a-th 
electronic bath and G^'^{iujn) that of the correspond- 
ing impurity problem, a local self-energy is extracted 
from the solution of the Anderson model for each sub- 
lattice a via EQ,,„(iw„) = G^ ,^(«w„)~i - G^^/iiuJn)"'^- 
The lattice Green's function within the dynamical mean- 
field approximation then reads Go-(k, iaj„) = [{iujn + 
/i)I — Eh — t(k) — Scr(iw„)]~^, where bold letters in- 
dicate N X N matrices in sublattice space, I is the 
identity matrix and t(k) is the Bloch matrix connect- 
ing the different sublattices. Matsubara frequencies, 
ujn — (2n -I- l)7r//3, with P = 200 are used. Specializing 
to the three- fold order realized on the triangular lattice, 
the matrix elements between sublattices a = A,B,C are 
tAsik) = tBc{k) = 4>k, tAcik) = (t>l, having defined 
(/)fc = ^j g_^^i-k ^j^fj Jj (i = 1^ 3) the vectors connecting an 
A site to its three B neighbors (see Fig. [T]-b). E// is a 
diagonal matrix whose elements EH,a = 3 '^^-^^ ^('^7) 
are, for each a, the mean-field electrostatic potentials 
of charges on the nearest-neighboring sites. The prob- 
lem is solved self-consistently by imposing the condition 
(1/^.) Ek Gaa,.(k, iw„) = Gtfiiujr^), with A^, the total 
number of lattice sites. We use Lanczos diagonalization 
as a problem solver, taking 10 auxiliary sites for each 
impurity problem. The spectral density for electrons on 
the a-sublattice is obtained as AQ_cr(w)=-Im Gqq^o- (w)/7r. 
We restrict our analysis to solutions where the B and C 
sublattices are equivalent. Practical calculations for a 
density n of holes are performed by taking n electrons 
per site and changing the sign of t in Eq. (1). 

The pinball liquid. Fig. [T]-a shows the phase diagram 
of the EHM obtained from DMFT on the triangular lat- 
tice at n = 1/2 (quarter filling), which basically agrees 




EHM, n^ 
PAM./u 




^^ 


"~^ 






d) 


PAM 


\ 


EHM 






6 


8 


10 


12 


1- 



u 



FIG. 1: (Color online) Pinball liquid phase on the triangular 
lattice from DMFT. (a) The (/-V phase diagram of the EHM 
on the triangular lattice at quarter filling, (b) Sketch of the 
pinball phase with pins on the A sublattice (red, large disks) 
and balls (blue, small) on the B — C sublattice. The shaded 
area is the unit cell, (c) and (d) are the charge occupations 
and quasiparticle weight Za as a function of U . Dashed lines 
are obtained from the PAM Eq. ((Sjl . 



with previous results obtained from exact diagonalization 
(ED) on small clusters lOllUl- In the weakly correlated 
limit, [/ — > 0, a transition from a homogeneous metal 
(HM) to a metallic CO state is realized upon increasing 
the intersite repulsion V ^ similar to what is found on bi- 
partite lattices [5ll7H9|. In contrast, the charge ordered 
insulator obtained on bipartite lattices in the large U re- 



gion is replaced here by a metallic CO, termed the pinball 
liquid 18|, ll9| , extending up to increasingly large values 
of V . In such phase, quasi-localized charges {pins) coex- 
ist with a low density of itinerant carriers (balls) flowing 
in the honeycomb lattice formed by interstitial sites (see 
Fig. [l}b). At large V and small U, the full CO in the 
three-fold geometry corresponds to a density ha — >■ 3/2 
on the charge-rich sites and ns — > on the charge-poor 
sites. Increasing U tends to suppress the double occu- 
pancy and eventually when U > 3^, the local Coulomb 
energy cost is suppressed by leaving only one charge per 
A site (a pin, ua — 1). The excess charge density is then 
transferred to the charge-poor sites (balls, ub — 1/4). 
The transition is sharp, as shown in Fig. [T]-c, and it 
can be anticipated that the two species behave very dif- 
ferently; pins form an almost half-fllled system with a 
marked Mott behavior due to strong onsite Coulomb in- 
teractions, in opposition to the balls which, owing to their 
low concentration, are in principle protected from them. 
We note here that the stability of the pinball phase is not 
limited to the case n = 1/2 described here, but actually 
extends to the entire range 1/3 < n < 2/3 20] . 

The quasiparticle weight of the correlated pins, defined 
as Za = [1 — dImT,A{^)/duj]~^Q (depicted in Fig. [Ijd) 
is large in the conventional CO phase at low U, and 
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FIG. 2: (Color online) Metallic behavior of the pinball liquid 
vs. square lattice charge ordered insulator. Spectral functions 
^^(a;) and ^s(cj) projected on the charge rich and charge 
poor sublattices, for [/ = 12 in the EHM on the triangular 
lattice at V^ = 3 (a) and on the square lattice sd, V = 2.5 in 
the COI phase (c). (b) Same quantities for the PAM at the 
corresponding value A = 4.5i. The vertical lines denote the 
Fermi energy. 



drops drastically upon entering the pinball liquid (here 
at C/ = 8.5t). As mentioned above, this occurs when 
the A-sublattice concentration spontaneously adjusts to 
riA — I (Fig. [B-c). Within the accuracy of our method, 
a small non- vanishing quasiparticle weight persists at all 
U, indicating that the pins do not undergo a full Mott 
transition, but remain in a strongly correlated Fermi liq- 
uid state up to f7 — > cx). We believe that this originates 
from the fact that the filling of the correlated hybrid 
band is never exactly integer. This contrasts with the 
situation encountered on bipartite lattices [Of , where the 
electron concentration in the correlated band is automat- 
ically fixed to half-filling once that the Fermi energy falls 
inside the gap, resulting in a true Mott transition for 
U>W a.ndV>2. 

The sublattice spectral densities ^0(0;) within the 
charge ordered pinball phase, illustrated in Fig. [2]-a for 
[/ = 12 and V^ = 3 for the EHM, confirm this pic- 
ture. The system is metallic, as testified by the strong 
spectral weight crossing the Fermi level, which primar- 
ily originates from the ball sublattice (B,C sites). The 
A-sublattice instead shows typical features of strongly 
correlated systems with a spectral weight mostly located 
in two separate regions away from the Fermi level - the 
lower and upper Hubbard bands separated by an energy 
U - and a weak quasiparticle residue at the Fermi energy. 



quantified by Za- The spectrum of the charge ordered 
insulator (COI) phase on the square lattice is shown in 
Fig. [2]-c for comparison. 

Periodic Anderson model and Emergent heavy 
fermions. To get further insight into the effects of 
electronic correlations we now map the problem into a 
periodic Anderson model (PAM). This is achieved by in- 
troducing creation operators af , bf and cf respectively 
on the A, B and C sublattices and diagonalizing the 
bands of the honeycomb lattice formed by B and C sites 
(tight-binding on graphene). Rewriting (j)k = J2i^^^''' 
as Eke^^'', we can define B-C subspace spinor opera- 
tors ^^j_ = (e^''^''/^,±e'^''/^)/-\/2, whose energies are 
ek,± — ^Ek (the origin is fixed at the Dirac point). 
Dropping the Coulomb interaction on the conduction 
bands leads to 



HpAM = ^ £Aal„ai„ + ^ ek,alkacr'^kacr (2) 

icr kaa 

+ ^{Vk,ae'^'''^'al^'jkaa + h.c.) + U'^al^a^-^al^aii 



where the hybridization between the conduction bands 
and the localized level is determined by Vk.^ = 
-V2EkSin{3ek/2) and Vk,+ = V2EkCOs{Wk/2). The 
energy of the localized level is set by the charge transfer 
gap, EA = —A, which originates from the electrostatic 
interaction between charges on different sublattices. The 
pinball phase studied in the preceding paragraphs cor- 
responds to A = 3V/2 and is therefore in the charge 
transfer regime U ^ A (because U > 3V ^^ 3V/2, cf. 
Fig.[T]-a). The quantities Ua, Za and ^0,(0;) calculated in 
the PAM are compared with the EHM results in Figs. [I]- 
(c,d) and[2]-(a,b) showing a consistent agreement in the 
pinball liquid phase. 

The nature of the coupling between pins and balls is 
now conveniently analyzed based on the low energy ef- 
fective PAM varying the number of electrons per unit 
cell ncoii G [0, 3] (riccii = in) for a fixed A — 4.5. This 
interval includes the region of stability of the pinball in 
the original EHM ^, i.e. 1 < riceii < 2. The density of 
states (DOS) of TJpam in the non-interacting limit U — 
is shown in Fig. [3]-a. Hybridization broadens the local- 
ized pin level at e^ = — A (vertical line in Fig.[3]-a) into a 
dispersive band, and also affects the dispersion of the con- 
duction bands, that is no longer symmetric around the 
Dirac point w = 0. Because ea < 0, pins are populated 
first and n^ — ?^ccii until ncoii = 1 as seen in tIq, depicted 
in Fig. [3]-b. For riceii > 1, the A sublattice is filled with 
~ 1 electron per site (the maximum allowed by the large 
U) and the remaining electrons move to the conduction 
band. Correspondingly, the Fermi energy moves from the 
bottom of the (— ) band to the Dirac point, reached at 
"-cell = 3. Interestingly, Za (blue squares in Fig. |3]-b) is 
strongly suppressed in the whole pinball phase (hatched 
region in Fig. |31-b). Following Refs. 25|, |26|, we define 
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FIG. 3; (Color online) Strong quasiparticle renormalization 
in pinball liquid phase and antiferromagnetic Kondo correla- 
tions, (a) DOS of the noninteracting PAM (!7 = 0) Q for 
A = 4.5. Different colors correspond to projection on the non- 
hybridized bands: pins (blue), — (orange) and + (red), (b) 
A/B sublattice densities (fuU lines=DMFT, dashed=ED) and 
pin quasiparticle weight (DMFT) as a function of the filling 
UcEii (c) hybridization, quasiparticle renormalization (DMFT) 
and magnetic correlations (ED). 



the effective self-energy for conduction electrons as 



^k,± 



{iUn) 



\Vk,: 



+ H~ EA- SA(iw„) 



\Vk3 



(3) 



iuj„+fi-Ek 



which accounts for hybridization as well as indirect inter- 
action effects arising via the Coulomb repulsion of local- 
ized charges. Based on Eq. ([3|). we evaluate the quasi- 
particle renormalization Zk,± = [1 — dlm{l]k,±) / dio]'^^Q 
of the balls. Note that Vk^- vanishes for wavevectors k 
forming an angle 7r/6 [tt/S] with the x-axis, leading to 
cold spots where conduction and localized electrons are 
perfectly decoupled and Z^ _ = 1 identically. In Fig. [3]- 
c we report instead Z^^- in the directions of maximum 
hybridization on the Fermi surface (hot spots), for which 
a sharp drop is found at nceii = 3/2. This is attributed 
to an increase in the hybridization with the correlated 
pins occurring around this filling (see also |T4^_p). On 
the same figure is also reported the inter-species mag- 
netic correlation (m^mf ) = ((n^t ~ 'T-yii)(^Bt ~ "-Si)) 
calculated by ED on a 12-site cluster (we use ED be- 
cause inter-site correlations cannot be addressed in the 
single site DMFT scheme). As shown in Fig. [3]-c, the 
strong renormalization of the metal is accompanied by 
a buildup of anti-ferromagnetic correlations, correspond- 
ing to a screening of the localized moments by the con- 



duction electrons. From our data we conclude that pins 
and balls are nearly decoupled for nceii < 3/2, while a 
strong Kondo coupling between the two fluids, respon- 
sible for strong renormalization of the conduction elec- 
trons, is present at larger concentrations. 

Wigner crystal in the 2DEG. According to the lattice 
analogy proposed in Refs. [y, [Zl, [ij], the existence of 
pronounced short range order close to the MIT enables a 
mapping of the continuous 2DEG to our model iJpAM at 
nceii = 1- We can then associate A oc r^ [ij], where rg is 
the ratio of Coulomb to kinetic energy in the electron gas. 
The case A = 4.5 analyzed in Fig. |3]-b then corresponds, 
for ncou = 1, to a 2DEG well into the Wigner crystal 
phase: the system is insulating, as both Za and the DOS 
at the Fermi eriergy (not shown) vanish at this point. 



Following Ref. [lj| the quantum melting of the Wigner 



crystal upon reducing Vg can be identified with a closing 
of the gap between the pin and (— ) bands (Fig. [3]-a) 
occurring below a critical Ac. From our DMFT results, 
this happens at Ac ~ 3.7 for U = f 2. Although a detailed 
study of the MIT of the 2DEG is beyond the scope of this 
work, we note that a metallic state can also be achieved 
at A > Ac, by doping away from riceii = 1- The resulting 
metal is very different when doping holes into the pin 
band (uccw = 1 + 6 with (5 < 0, leading to Za oc |(5|, cf. 
Fig. [3]-b) or electrons into the conduction band {5 > 0, 
with Za jumping sharply to a small but finite value). The 
latter situation actually corresponds to the self-doping 
instability proposed in Ref. [6|. If this scenario were 
realized, the metal on the verge of Wigner crystallization 
would actually be a pinball liquid itself, indicating an 
appealing connection between the physics of the 2DEG 
and the emergence of heavy fermion behavior. 

Concluding remarks. The combination of electronic 
correlations and geometrical frustration on charge or- 
dered triangular lattices leads to an exotic pinball liq- 
uid state in which localized and itinerant electrons co- 
exist. Based on the close analogies with heavy fermion 
systems, pointed out via an explicit mapping of the elec- 
tronic problem onto the periodic Anderson model, non- 
Fermi liquid behavior should occur in such phase as a 
consequence of the scattering of the itinerant electrons off 
the localized moments of the pins. On the experimental 
side, this behavior would be consistent with the anoma- 
lous metallic phases of the family of organic conductors 
e'-(BEDT-TTF)2X which display violation of the Mott- 
loffe-Regel limit in the resistivity, absence of Drude peak 
in the optical conductivity 22l 12311 and non-monotonic 
1/TiT-NMR relaxation rates [28j close to charge or- 
dering. The senario proposed here for emergent heavy 
fermion behavior could also be of relevance to other cor- 
related materials exhibiting charge ordering on trian- 
gular lattices such as transition-metal dichalchogenides 
EEEI. AgNi02 [29] and ^He bilayers [30|, as well as 
to the interaction-driven MIT in the 2DEG 91 . 
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